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In this work, we present a new model for the heat conductivity of porous dust layers in vacuum, based on an existing solution 
of the heat transfer equation of single spheres in contact. This model is capable of distinguishing between two different types of 
dust layers: dust layers composed of single particles (simple model) and dust layers consisting of individual aggregates (complex 
model). Additionally, we describe laboratory experiments, which were used to measure the heat conductivity of porous dust layers, 
in order to test the model. We found that the model predictions are in an excellent agreement with the experimental results, if we 
include radiative heat transport in the model. This implies that radiation plays an important role for the heat transport in porous 
materials. Furthermore, the influence of this new model on the Hertz factor are demonstrated and the implications of this new 
model on the modehng of cometary activity are discussed. Finally, the limitations of this new model are critically reviewed. 
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1. Introduction 

Small bodies in the Solar System are either covered by re- 
Q golith, a granular medium formed through the continuous bom- 
bardment by interplanetary bodies (asteroids) or even consist 
entirely of granular material (comet nuclei). The heat conduc- 
ed tivity of the surface layers of such bodies plays an important 
role for the minimum and maximum temperature as well as 
CT^ for the temporal temperature evolution on the surfaces of these 



> 

in 
cn 
in 
o 



bodies. It also determines the amount of energy transported into 
deeper regions of these bodies and, thus, the temperature strat- 
ification. In the young Solar System, the fate of planetesimals 
with respect to sintering and melting is also determined by the 
heat conductivity of the granular material they consist of (see, 
e. g., |Hevey and Sanders} |2006[ [Gupta and Sahijpal} |2007 



[Moskovitz and Gaidos 2011[ " Krause et al. 201 1 1. For comet 
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nuclei, the heat conductivity of the ice-free surface layers de- 
termines (among other factors) the temperature and, thus, the 
evaporation rate of the underlying volatiles. 

In paper I ( |Gundlach et al. 2011bj , we presented measure- 



ments of the gas transport through porous dust layers and de- 
rived the gas permeability for arbitrary layers of (spherical and 
monodisperse) SiOa dust particles. In this paper (paper II), 
we describe laboratory experiments on the determination of the 
heat conductivity of granular dust layers consisting of spherical 
Si02 dust particles of different sizes (Sect. |2]-|4|i. Thereafter 
(Sect. [5]), we adopt a model for the heat conductivity of packed 
spherical granular particles by accounting for the mutual van 
der Waals interaction of the grains and extend the model to ar- 
bitrary granular packing of dust aggregates, as hypothesized to 
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exist on the surfaces of comet nuclei ( |Skorov and Bluml|201 l| l. 
In Sect. |6j we discuss the implications of this model on the 
Hertz factor and on the modeling of cometary activity. Further- 
more, the limitations of this new model are critically reviewed. 
Finally, we conclude this paper by critically addressing the ca- 
pabilities of the heat-conductivity model for future applications 
in planetary science (Sect. |7]i. 

2. Experimental technique 

2.1. Setup 

The experiments were performed using the experimental 
setup introduced in paper I (see Fig. [TJ. Minor modifications of 
the setup were conducted in order to measure the surface tem- 
perature and, therewith, the heat conductivity of porous dust 
layers. 

The measurements were performed in a high-vacuum cham- 
ber (see no. 1 in Fig. [T]i at pressures below 10"'*mbar. In this 
case, the mean free path of the gas molecules is much larger 
than the mean pore size of the dust layer Thus, heat conduc- 
tion due to gas diffusion can be neglected in the experiments. 

The dust samples (2; see Sect. |3]l were filled into a cylindrical 
sample holder (3), which was positioned on the cold plate (4). 
In paper I, the cold plate was cooled by the cooling system (5), 
which was not used in this work. Here, the cold plate acts as a 
constant-temperature reservoir 

The dust samples were irradiated with a halogen lamp (6; in- 
tensity: 5.32 solar constants on the surface of the sample). The 
beam profile of the light source was adjusted with two lenses 
(7). Different grey filters (8) with transmission rates of 50 %, 
20 %, 10 % and 5 %, were introduced into the light beam to vary 
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1: High-Vacuum Chamber 

2: Dust Sample 

3: Sample Holder 

4: Cold Plate 

5: Cooling System 



6: Halogen Lamp 
7: Lenses 
8: Grey Filter 
9: Alumium Plate 
10: Infrared Camera 



11: Germanium Window 
12: Temperature Sensor 
13: Computer System 



Figure 1 : Schematic diagram of the adapted experimental setup. See |Gundlach| 
|et al.|j2011b) for the unmodified version. 



the intensity of the irradiation, without changing the spectrum 
of the Hght source. 

An aluminum ring (9) was implemented on top of the sample 
holder, to avoid additional heating of the sample holder and, 
therewith, additional heating of the dust samples. In order to 
minimize the heat flow through the sample holder, a material 
with a low specific heat conductivity (for comparison, metals, 
glass or plastic materials typically have specific heat conductiv- 
ities of /I = 10 - 400WK-im-i, A = 1 - SWK-'m-i and 
A - 0.2 - 0.5WK 'm~') was chosen (polyvinyl chloride, 
A - 0.15 WK' m ') and the thickness of the wall of the sam- 
ple holder was small (2 mm) in comparison with the diameter 
of the dust samples (25 mm) so that the amount of heat trans- 
ported through the sample holder is much smaller than the heat 
transported through the dust sample. Furthermore, the sample 
holder was in a very loose contact with the cold plate, which 
further decreased the heat flow through the sample holder. 

The surface temperature of the dust layer was measured with 
an infrared camera (10) through a germanium window (11). 
Due to the low transmissivity of germanium in the infrared 
(transmission coefficient; ~ 0.4), a calibration of the infrared 
camera observing the surface of the dust sample through the 
germanium window was necessary (see Sect. 2.2 1. A temper- 



ature sensor (12) was incorporated into the cold plate, directly 
beneath the dust sample (distance between the bottom of the 
dust sample and the temperature sensor: 2 mm). With this sen- 
sor, the temperature evolution of the cold plate and, therefore, 
the temperature evolution of the bottom of the dust sample was 
monitored. The temperature of the cold plate remained con- 



stant during the experiments, due to the high heat capacity of 
the material (see Fig. |4]in Sect. |4|l. To control the experiment, 
a computer system (13) was used. 

2.2. Calibration of the infrared camera 

Observing the surface of the dust sample with the infrared 
camera through a germanium window leads to a falsification of 
the observed temperature. Thus, a calibration of the infrared 
camera observing the surface of the sample through the germa- 
nium window was necessary. Therefore, the temperature evolu- 
tion of a stove top was measured with and without germanium 
window. From this calibration measurement, the influence of 
the germanium window on the measured temperature was de- 
rived. The difference between the measured temperature with 
and without germanium window was used to calculate the real 
surface temperature of the dust sample. 

2.3. Procedure 

At the beginning of the experiments, the dust samples were 
filled into the sample holder on top of the cold plate. Then, 
the experimental chamber was evacuated to low pressures (p < 
10""* mbar), so that the mean free path of the gas molecules was 
much larger than the pore size of the dust sample. In this case, 
heat transport due to the gas diffusion was negligible. After that, 
the dust samples were irradiated by the light source until the 
equilibrium temperature was reached (see Fig. |4]in see Sect. |4]l. 
The temperature evolution of the surface of the dust sample and 
the bottom of the dust sample were measured by the infrared 
camera and the temperature sensor, respectively. 

3. The dust samples 

Five different dust samples (DSl, DS2, DS3, DS4 and DS5) 
were used in the experiments. All samples consisted of spher- 
ical, black Si02 particles (see Fig. |2|l. In this section, we dis- 
cuss the physical properties of the dust samples relevant for this 
work, such as the size distributions of the particles, the vol- 
ume filling factors of the dust samples, the plane albedo and the 
emissivity of the dust samples. 

3.1. Size distribution and volume filling factor 

For the estimation of the particle diameters s, we used a light 
microscope (DSl, DS2, DS3 and DS4) and a cafliper (DS5). 
Fig. [3] shows the resulting particle diameters. The volume fill- 
ing factors 4> of the dust samples were measured by comparing 
the mass of the dust sample with the mass of solid SiOa, which 
would occupy the same volume. 

Tab. [T] summarizes the mean particle diameters and the vol- 
ume filling factors of the five different dust samples. 

3.2. Albedo 

The plane albedo of the dust samples was estimated using 



the experimental technique invented by Krause et al. (201 1 1. In 
this experiment, porous dust samples were positioned inside a 
vacuum chamber and irradiated by an infrared-laser. The spa- 
tial and the temporal temperature evolutions of the surface of 
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Figure 2: Images of a black dust sample (a; DSl) and of the individual SiOi particles: DSl (b), DS2 (c), DS3 (d), DS4 (e) and DS5 (f). The images a and f were 
taken with a digital camera and the images b - e were taken using a light microscope. 



Table 1: Comparison of the mean particle diameters s and the volume filling 
factors (p of the five dust samples (DSl, DS2, DS3, DS4 and DS5). The errors 
of the mean particle diameters and the volume filling factors are given by the 
standard deviations of the measurements. 



Dust Sample 


s Lum] 






DSl 


40 ± 15 


0.670 ± 0.015 


DS2 


177 ± 59 


0.634 ± 0.024 


DS3 


371 ± 126 


0.68 


8 ±0.018 


DS4 


507 ± 84 


0.685 ± 0.017 


DS5 


925 ± 64 


0.68 


8 ± 0.008 



the dust sample were monitored with an infrared camera. The 
heat conductivity and the plane albedo of the dust samples were 
independently derived by comparing the experimental results 
with a thermophysical model. 

Using this technique, a plane albedo of A - 0.35 was esti- 
mated for the dust samples (Krause and Skorov, personal com- 
munication). 

3.3. Emissivity 

The determination of the emissivity of the dust sample is im- 
portant for the calculation of the heat conductivity using the 
energy balance equation. Thus, the emissivity of the dust sam- 
ples was measured with the infrared camera. The samples were 




10 100 1000 

Particle Diameter [t<m] 



Figure 3: Particle size distributions for the five difterent dust samples (solid 
curves). The mean particle diameters are visualized by the dashed lines: i = 
(40 ± 15) /<m (DSl), s = {111 ± 59) /^m (DS2), s = (370 ± 126) yum (DS3), 
s = (507 ± 84) /im (DS4) and s = (925 ± 64))jum (DS5). 

not heated or cooled so that the temperature of the sample was 
equal to the ambient temperature, which was monitored by sev- 
eral temperature sensors. For this measurement, the observation 
angle between the surface of the sample and the infrared cam- 
era was 45°, which is identical to the observation angle in the 
experiments. 

The software of the infrared camera was then used to ad- 
just the emissivity until the infrared camera has measured the 



3 



T" 



"T" 



500 



Surface Temperature of the Dust Layer 
Surface Temperature of the Sample Holder 
Bottom Temperature of the Dust Layer 




200 



400 600 
Time [s] 



800 



^ 550 

Z. 500 

tn 

Q 
0) 

£ 450 



o 


Q. 



ID 
U 



400 



350 



300 





h 


^ 5.74 mm 


A 


h 


= 4.38 mm 


X 


h 


= 2.36 mm 


+ 


h 


= 1 .30 mm 


□ 


h 


= 0.74 mm 



It 



X 



I 
X 

...a- 



2 3 4 

Intensity / Solar Constant 



Figure 4: Temporal evolution of the surface temperature of the dust layer (solid 
curve), of the surface temperature of the sample holder (dotted curve) and of the 
bottom temperature of the dust layer (dashed curve) during and after iiTadiation 
(/ = 5.32/s). In this experiment, the DSl dust sample with a height of ft = 
(1.30 ± 0.01)mm was used. The gray regions along the curves represent the 
measurement uncertainties. 



Figure 5: Equilibrium surface temperatures of the dust layers as a function of 
the intensity of the light source. The experiments were carried out for five dif- 
ferent dust layer heights, (0.74 ±0.01) mm (squares), (1.30 ±0.01) mm (pluses), 
(2.36 ± 0.01) mm (crosses), (4.38 ± 0.01) mm (triangles) and (5.74 ± 0.01) mm 
(asterisks). The dotted lines are introduced to guide the eye. 



correct temperature of the dust sample, assuming Lambertian 
emission. Using this technique, an emissivity of soust — 0.99 
was derived. Here, the emissivity and the plane albedo (over 
all wavelengths) can not be correlated via Kirchhoff's law, be- 
cause the emissivity was only measured within the spectral 
range of the infrared camera, limited by the germanium win- 
dow (7.5 jum < wavelength < ~ 12 /im). 



4. Experimental results 

4.1. Surface Temperature 

We measured the temporal evolution of the surface temper- 
ature of the dust layer, Tsurf, and of the bottom temperature 
of the dust layer, T Bottom, for different heights of the dust layer, 
h, different particle sizes and different intensities of the light 
source. The real surface temperature was calculated using the 
calibration of the infrared camera, which was discussed in Sect. 
OandllJ] 

Fig. |4]shows the temporal temperature evolutions of the sur- 
face of the dust layer (solid curve), of the surface of the sample 
holder (dotted curve) and of the bottom of the dust layer (dashed 
curve) during and after irradiation. In this example, the DSl 
dust sample with a height of h - (1.30 + 0.01) mm was used. 
The surface was irradiated with an intensity of / = 5.32 Ig (Ig : 
Solar constant). During all experiments, the bottom tempera- 
ture of the dust layer was Tsottom - ~ 300 K. Note that the 
change of the bottom temperature, NT Bottom < 2K, is small 
compared to the change of the surface temperature during the 
experiment. The error on the temperature measurements was 
tsTsui-f ± 2.0 K (infrared camera) and ArBo„o„, + 1.5 K (tem- 
perature sensor). The errors are visualized by the gray areas in 
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Figure 6: Equilibrium surface temperatures of the dust layers with respect to 
the height of the dust layer for dift'erent intensities of the light source: 0.27 1$ 
(squares), 0.53 /s (pluses), 1.06 (crosses), 2.66 /s (triangles) and 5.32 (as- 
terisks). The dotted lines are introduced to guide the eye. 



In total, 25 experiments were carried out using the DS 1 dust 
samples. Five different dust layer heights, (0.74 ± 0.01) mm, 
(1.30 + 0.01) mm, (2.36 + 0.01) mm, (4.38 + 0.01) mm and 
(5.74 ± 0.01) mm, and five different intensities of the light 
source, 0.27 Is , 0.53 Is , 1 .06 Is , 2.66 Is and 5.32 Is , were used. 
Fig. |5] and |6] are showing the resulting equilibrium tempera- 
tures, obtained in these experiments, with respect to the inten- 
sity of the light source and with respect to the height of the dust 
layer. The dotted lines are implemented to guide the eye. 

Furthermore, the evolution of the surface temperature of the 
dust layer was measured for different particles sizes. Fig |7] 
shows that the equilibrium surface temperature of the dust layer 
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Figure 7: Resulting equilibrium surface temperatures of the dust layers for dif- 
ferent mean particle diameters. In these experiments, the height of the dust lay- 
ers and the intensity were kept constant, A = (2.36±0.01)mmand/ = I.O6/5, 
respectively. 

increases slightly with the mean diameter of the particles. In 
these experiments, the height of the dust layers and the inten- 
sity were kept constant, /; - (2.36 + 0.01) mm and / = 1.06 Is- 
The mean diameters of the particles were summarized in Sect. 

4.2. Energy balance 

The energy balance equation was used to calculate the heat 
conduction through the dust layers in thermodynamic equilib- 
rium (constant surface temperature), when the different energy 
fluxes remained constant. Then, the conducted heat though the 
dust layer is given by, 

Econd - Ejrr - Egad , (1) 

where Ejrr - - A) is the energy input per unit surface 
into the sample due to irradiation and Enad = cr ( eoust T^urf ~ 
SBack T'^aci) thermal radiation of the surface of the dust 

layer, including the thermal radiation of the background. I is 
the intensity of the light source, A is the plane albedo of the 
dust sample, cr is Boltzmann's constant, eoust is the emissivity 
of the dust sample, e^ack is the emissivity of the background 
and Tsack is the temperature of the background, with TBack - 
300 + 2 K (ambient temperature). In this work we assume that 
the emissivity of the background is unity (black body radiation). 

We calculated the conducted heat through the dust layer for 
the 25 experiments, in which the DS 1 dust samples were used. 
Fig. |8]shows the conducted heat, normalized to the energy input 
into the sample, with respect to the height of the dust layer. The 
normalized conducted heat decreases with the height of the dust 
layer, which can easily be explained with the decrease of the 
temperature gradient inside the layer. The temperature gradient 
decreases linearly with the height of the dust layer, if the sur- 
face temperature and the bottom temperature of the dust layer 
are fixed. If the surface temperature increases with increasing 



Figure 8: Conducted heat through the dust layers, normalized to the energy 
input into the samples. The heat fluxes were calculated for the 25 experiments, 
in which the DSl dust samples were used (asterisks, triangles, crosses, pluses 
and squares). A good fit to the data (solid curve) is given by Eq. [2]with a = 
(7.41 ±0.29) mm-'. 

height of the dust layer, the decrease of the temperature gradient 
is smaller compared with the case discussed above. This behav- 
ior can be found in Fig. |8] because the maximum surface tem- 
peratures of the dust layers were still slightly increasing with 
increasing heights of the dust layers (see Fig. [6]l. However, 
a transition to the linear decrease of the normalized conducted 
heat through the dust layer can also be seen (see Fig. |8]l. Thus, a 
reasonable fit to the data (solid curve) is given by the following 
function, 

Eirr 1 + fl"' h ' 

witiia = (1.35 + 0.05)mm. 

The error due to heat conduction through the sample holder 
can be neglected, because of the shallow temperature increase 
of the surface of the sample holder (see Fig. |4]) and because of 
the small wall thickness of the sample holder compared with the 
diameter of the dust layer Thus, the uncertainty of the derived 
heat flux trough the dust layer is mostly given by the uncertainty 
of the temperature measurement. 

4.3. Heat conductivity 

In porous materials, heat can be transported inside the solid 
material, due to radiation inside the pores of the material and 
due to gas diffusion (see Fig. |9]l. The total heat conductivity 
of porous materials is a then given by the combination of these 
three effects. 

Air, T, A(r)) = AsoUdir, T) + Ar,,{T, A(r)) + Acas ■ (3) 

Here, Asoiidir, T) is the heat conductivity through the solid net- 
work of particles and r is the radius of the particles. The heat 
conductivity of a porous material can be orders of magnitudes 
lower than the heat conductivity of the bulk material, due to the 
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Figure 9: Sketch of the heat transport in a porous dust layer. Energy can be 
transported inside the solid material (solid arrow), due to radiation inside the 
pores of the material (dashed an'ow) and due to gas diffusion (dotted aiTow). 
In this example, all energy fluxes are directed downwards, because the surface 
temperature of the dust layer is higher than the temperature of the bottom. 



small contact areas between the particles, which are determin- 
ing the efficiency of the heat exchange between the particles. 
This effect is often taken into account by the introduction of the 
so-called Hertz factor H(r, T), Asoiidir, T) = H(r, T) AB„ik(T) 
(see Sect. |6.1| for details), where AsuikiT) is the heat conduc- 
tivity of the bulk material. For vitreous Si02, the temperature 
dependence of the bulk heat conductivity can be described with 
a linear function. 



,(r) ^ biT + b2.. 



withal = (1.26 ± 0.05)X IQ-^WK-^m' and 



0.20) X 10 ' WK' m-' (Ratcliffe 1963 1. 



(4) 
(9.94 ± 



The radiative heat transport in a porous material can be cal- 



culated treating the radiation as a photon gas (Merrill 1969 1, 



or by the formulation introduced by |Schotte| ( 19601. Thus, the 
heat conductivity of the radiative transport can be written as. 



with 



AR,AT,Mr)) ^ KT'A(r) (5) 

("16/3 0- (Merrill, 1969) 
I 8 ere (Schotte, 1960) ' 



where A(r) is the mean free path of the photons. Note that the 
first of the two formulations of the radiative heat conductivity 
does not depend on the emissivity of the dust layer In this work, 
we used Eq. |5]with k = 16 / 3 cr. In our case (eoust = 0.99, see 
Sect. |3.3| l, both formulations of the radiative heat conductiv- 
ity differ only by a factor 2/3, which we regard as negligible. 
However, for the calculation of the radiative heat conductivity 
for materials with an emissivity strongly deviating from unity, 
the second formulation of the radiative heat conductivity should 
be used. 



Furthermore, we can neglect heat transport due to gas diffu- 
sion A(jas = 0, because of the low pressures during the experi- 
ments (see Sect. |2.3| l. 

We calculated the heat conductivity of the dust samples as- 
suming a constant vertical temperature gradient inside the dust 
layers in the thermodynamic equilibrium. This assumption only 
holds if the heat conductivity is not dominated by the radiative 



heat transport (see Sect. 6.2 1. Here, we used relatively compact 
samples (0 - 0.63 - 0.69, see Sect. |3.11l, where the radiative 



heat transport is present (see Sect. [5]T i, but not dominating the 
heat transport. In this case, the heat conductivity can be calcu- 
lated as follows. 



A{r,T,A(r)) 



-'Cond' 



(h) 



'Surf 



- T„ 



(6) 



Fig. 10 shows the derived heat conductivities of the DSl dust 



samples (see also Fig. |5]and Fig. |6]l with respect to the surface 
temperature of the dust layers. The heat conductivity increases 
with the surface temperature of the dust layer The uncertainty 
of the calculated heat conductivities is primarily given by the 
error of the temperature measurement. Due to the usage of Eq. 
[T] the derived heat conductivity is calculated from the difference 
of two large numbers. As a result, the scatter of the derived 
heat conductivities is bigger for lower temperatures {Tsurf ^ 
Tsack) than for higher temperatures (Ts„rf > TBack)- For this 
reason, we have not included the experiments performed with 
grey filters having very low transmission rates of 1 % and 2 %. 
In Sect. 5.1 the derived heat conductivities are compared with 



a general model for the heat conductivity of porous dust layers. 

Furthermore, the influence of the grain size on the surface 
temperature of the dust layers and, therewith, on the heat con- 
ductivity was investigated (see Fig. |7]). Therefore, the height of 
the dust layers was kept constant, h - (2.36 ± 0.01) mm. Fig. 
1 1 shows the derived heat conductivities (squares) with respect 
to the mean particle diameters. The heat conductivity decreases 
with increasing mean particle diameter. Note that the heat con- 
ductivity is a combination of the heat conductivity of the solid 
material and the heat conductivity of the radiative transport. 

Additionally, the heat conductivities of porous dust samples 
composed of spherical, 1.5 micrometer-sized SiOa particles 
( [Krause et al.| |201 l| l were implemented in Fig. [TT] (crosses). 
The lower value, A - 0.021 WK ' m denotes the heat con- 
ductivity of the dust layer with the highest volume filling factor. 



(/i = 0.54, measured by Krause et al. (201 1 1. The upper value, 
A = 0.043 WK^' m^', was derived for = 0.67, using the 
correlation between the heat conductivity and the volume fill- 
ing factor of the material ( Krause et al.||20lT| l, 



Aicfi) - ci exp [ C2 < 



I , (7) 
1 andc2 = 5.19±0.73 



wiflici = (1.31±0.48)xlO-^WK-'m 
(Krause, personal communication). These fit parameters were 
obtained by fitting the data obtained by Krause et al. ( 201 l| l 
without using the value for fused silica. 

A detailed theoretical discussion of the heat conductivity of 
porous dust layers is given in Sect. |5] where a general model for 
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Figure 10: Temperature dependence of the heat conductivity (DSl dust sam- 
ples; diameter: (40 ± 15)jum), measured for different dust layer heights: 
(0.74 ± O.OI)mm (squares), (1.30 ± 0.01) mm (pluses), (2.36 ± 0.01) mm 
(crosses), (4.38 ± 0.01) mm (triangles) and (5.74 ± 0.01) min (asterisks). The 
dashed curve shows the temperature dependence of the heat conductivity cal- 
culated with the modified model of ]Chan and TienH1973| see Sect.[5l]Eq.[T2) 
assuming that the specific surface energy of Si02 is not temperature dependent, 
ysiOi = 0.02 J m"~. The model predictions for a temperature dependent spe- 
cific surface energy (see Eq.[TTJ are visualized by the solid curve. Additionally, 
the influence of the radiative heat transport (see Eq. |5j is shown for different 
values of the mean free path, 0.5 r, 1.0 r, 1.5 r, 2.0 r and 2.5 r, by the dotted 
curves. For the theoretical derivations of the heat conductivities, the tempera- 
ture inside the dust layer was averaged, T = T Bottom + ( Tsurf - T Bottom) 1 2 
(see Sect. |4.3|for explanation) 



the conductivity of dust layers will be presented and compared 
with the experimental results. 

5. A model for the heat conductivity in loose granular me- 
dia 

5.1. A general model for the heat conductivity of dust layers 

A very detailed analytical study of the heat conductivity of 
packed spheres in vacuum was carried out by Chan and Tien| 
( 1973| l. They solved the heat transfer equation of single spheres 
and used this result to calculate the conductivity of packed 
spheres. In their model, heat can be transported from one sphere 
to an other only through the contact interface, which can be de- 



termined using the Hertzian theory of contact (Hertz 1881 
They found that the heat conductivity of the packed spheres in 
vacuum can be calculated as follows. 



with 



^Bulk 



(T) 



3 I - 

4 ^(r) 



F(r)r 



1/3 



1 



0.5315(0) Niir) 



(8) 



(9) 



Here, f^ and £(7") are the Poisson's ratio and Young's modu- 
lus of the material, respectively. 5(0) is a model parameter. 



Figure 1 1 : Heat conductivity as a function of the mean particle diameter The 
experimental results are shown by the squares. In all experiments, the height of 
the dust layer was identical, h = (2.36 ± 0.01) mm. For comparison, the results 
of the model for the heat conductivity of packed spheres with </> = 0.68 in vac- 
uum (see Eq. |12) are visualized without the influence of radiative heat transport 
(dashed curve) and with the influence of radiative heat transport (dotted curves; 
see Eq. [Sj. The calculations were performed for a mean temperature of the dust 
layer of 340 K and for five dilferent values of the mean free path of the photons 
inside the pores of the porous material, 0.5 r, 1.0 r, 1.5 r, 2.0 r and 2.5 r. The 
crosses are experiinental results from |Krause et al.]j2011^ for dust layers with 
packing densities of = 0.54 and = 0.67. 



which depends on the packing structure of the material. NAir) 
and Niir) are the number of particles per unit area and unit 
length, respectively. Tab. [2] summarizes the constants S (0), 
NA{r) and Niir) for three different compact packing structures, 
simple cubic (sc), body-centered cubic (bcc) and face-centered 
cubic (fee). 

F(r) describes the external force, which acts on the spheres 
and determines the contact interface between the particles. 
Chan and Tien ( 1973| l used an external applied load or the 
weight of the spheres above the contact to calculate the contact 
interface. However, the adhesive bonding between two spheri- 
cal particles in contact (van der Waals bonding) can be orders 
of magnitude higher than their weight, when no external load is 
applied. The adhesive bonding between two spherical particles 



Table 2: Summarization of the constants S, NA(r), Nt{r) and the volume fill- 
ing factors (j) of the thi'ce different packing structures, simple cubic (sc), body- 
centered cubic (bcc) and face-centered cubic (fee). See |Chan and Tien]jl973) 
for details. 



Parameter 


SC 


BCC 


FCC 


SW 


I 


1 

4 


1 

3 


Niir) 


1 

4, -2 
1 


3 

16/- 
3I/2 
2r 


2(3)1/2^2 
(3/8)1/2 





0.52 


0.68 


0.74 
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in contact can be calculated, using the theory of Johnson et al. 
( fT97T) , 

F,Mr,T) ^ 3ny{T)r, (10) 

where, y{T) is the specific surface energy of the material, which 
is a measure for the adhesive bonding strength of the material. 
For SiOa, measurements of the specific surface energy at room 
temperature yield a mean value of jsiOi - 0.02 J m"^ (|Kendall 



letaL] [T9871 [Heim et al.[ [T9991 IGundlachlTaLi [20TTal l. How- 
ever, the specific surface energy should also depend on temper- 
ature, because an increase of the temperature increases the rate 
of induced dipoles per molecule. Since the influence of the tem- 
perature on the specific surface energy is not known in detail, 
we assume that the specific surface energy increases linearly 
with temperature, 

ysio,(T) ^ diT , (11) 

with d\ = 6.67 x 10"^ Jm^^K"^. The slope is given by the 
assumptions that the specific surface energy should be reduced 
to zero if the temperature tends to zero and that the specific 
surface energy is 75,0, - 0.02 J m"^ for T = 300 K. 

In this work, we use the theory of 'Chan and Tien ( 1973 see 
Eq. |8]l to calculate the heat conductivity of our dust samples, 
but instead of taking the weight of the particles into account, we 
use the adhesive bonding force (see Eq. [TOjl, to determine the 
contact area between the particles. Thus, the heat conductivity 
of porous media in vacuum without external load and radiative 
transport is given by. 



E(T) 



ny{T) r 



^{r,(P). 



(12) 



The dust samples used in this work are random close packed, 
(f) - 0.63 - 0.69 (see Sect. 3.1 1. Since the coordination number 



of a random close packed material is very close to the coordina- 



tion number of a body-centered cubic packed material ( van de 



[Lagemaat et al. 2001 van Antwerpen et al. 2010| l, we use the 
body-centered cubic packing structure to calculate the heat con- 
ductivity of the dust samples (see Tab. [2|. 

The bulk heat conductivity and the specific surface energy of 
SiOa are given by Eq. |4]and Eq. 1 1 respectively. For complete- 
ness it should be mentioned that Young's modulus of SiOa is 
also temperature dependent (Marx and Sivertsen 1953 ';'SpinnerJ 
[and Cleek||1960| l. However, the variation of Young's modulus 
due to temperature changes is less than 10% for temperatures 
between 300 K and 600 K, which was the highest surface tem- 
perature achieved during the experiments. Therefore, we use 
a constant value for the Young's modulus and Pois son's ratio 
of Si02, Es io^ = 5.51 X 10'°Pa and//5,o, = 0.17 ( jChan aiid 
|Tienl[T973l l. 

Using Eq, 



12 the heat conductivity of packed Si02 spheres 



for different temperatures was calculated and compared with 
the experimental results (see Fig. 



lOi 



The heat conductiv- 
ity was derived assuming a constant (i. e. temperature inde- 
pendent) specific surface energy, ysio^ - 0.02 J m"^ (dashed 
curve). The model is in a good agreement with the experimen- 
tal results for temperatures below ~ 400 K. But, for higher 
temperatures, the model predictions deviate from the data. 
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Figure 12: Theoretical derivation of the mean free path of the photons inside the 
voids of porous dust layers for different volume filling factors of the material. 
The data (crosses; [Skorov et al.||20lT} were fitted with Eq. |13| (dotted curve) 
to extrapolate the mean free path to higher volume filling factors. 



Including a temperature dependent specific surface energy 



(see Eq. 1 1 1 yields a better agreement between the model (solid 
curve) and the experimental results {x^ decreases from 9.19 X 
10"^ to 6.43 X 10"^). Additionally, the influence of the radiative 
heat transfer was included (see Eq. |5]l, using diff'erent values 
for the mean free path of the photons inside the voids of the 
porous material, A(r) = 0.5 r (x^ - 6.39 x 10"''), 1.0 r ix^ = 
6.90x10-''), 1.5 rOt-^ = 7.94 x 10"''), 2.0 r = 9.51x10"'') 
and 2.5 r (x^ - 1.16 x 10"^; dotted curves). For A(r) = 0.5 r 
(X^ = 6.39 X 10"*), the model predictions are in best agreement 
with the experimental results. 

The mean free path inside the voids of a porous material 



was recently investigated theoretically by Skorov et al. (2011 



crosses in Fig. \V2\ . They found that the mean free path strongly 
depends on the volume filling factor of the material. Unfor- 
tunately, the highest volume filling factor used in this work, 
(p - 0.35, is two times lower than the volume filling factor of 
our dust layers (see Sect. 3.1 1. Thus, we used the following fit 
function, which was theoretical predicted by Dullien ( 1991 1, to 



extrapolate the mean free path derived by Skorov et al. ( 201 I| l 
to higher volume filling factors (see dotted curve in Fig. |12| l. 



A(r, 



1 - <p 



(13) 



The best fit to the data is given by ei = 1.34 + 0.01. Using 
this formula, we can estimate the mean free path of the pho- 
tons inside a material with arbitrary volume filling factor. For 
(p - 0.67, Eq. [T3]yields A(r) = 0.7 r, which is in a good agree- 
ment with the mean free path derived from the comparison of 
our experimental results with the model for the heat conductiv- 
ity. Note that our experiments were performed with polydis- 
perse Si02 particles (see Fig. |2]and[3]) while the calculations 
were conducted for monodisperse particles. The presence of 
different particle sizes should influence the mean free path of 
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the photons inside the pores of the material. Another uncer- 
tainty is given by the assumption that the specific surface en- 
ergy increases Hnearly with the mean temperature of the dust 
layer (see Eq. 
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Furthermore, Eq. 



121 was used to calculate the influence of 

n 



the grain size on the heat conductivity of the dust layer Fig 
shows a comparison of the derived heat conductivities (dashed 
curve) with the experimental results (crosses and squares). Note 
that the heat conductivity of dust samples composed of spheri- 
cal, monodisperse 1 .5 micrometer-sized particles (crosses) was 
measured by [Krause et al.| ( |201 l| l. See Sect. |4.3| for a more 
detailed description. 

The calculations were conducted using a mean temperature 
of 340 K for the dust layer. This value was used because the 
maximum surface temperature of the dust layers was ~ 380 K 
(see Fig. [Tji and the bottom temperature was ~ 300 K during all 
experiments (see Fig. |4|. However, this value can slightly de- 
viate from the real mean temperature of the dust layer, because 
of the temperature dependence of the heat conductivity, which 
leads to a non linear temperature gradient inside the dust layer 
(see Sect. |6.2| i. 

The influence of the radiative heat transport (see Eq. |5]) on 
the model for the heat conductivity is also shown in Fig. [TT] 
(dotted curves) for different values for the mean free path of the 
photons inside the voids of the porous material, A(r) - 0.5 r, 
1.0 r, 1.5 r, 2.0 r and 2.5 r (dotted curves). A good agreement 
between the experimental results and the model can be found 
for A(r) - 0.5 r, which is consistent with the theoretical pre- 
diction of Eq. 13 A(r) - 0.7 r. Note that the heat conductivity 



of the particle network decreases with increasing particle diam- 
eter. However, for large particle diameters, the radiative heat 
transport becomes important, due to the increase of the mean 
free path of the photons (see Eq. [T3| ). As a result, the heat 
conductivity increases with increasing particle diameter The 
combination of both behaviors leads to a minimum of the heat 
conductivity for a certain particle diameter, which depends on 
the temperature on the material. 

5.2. Heat conductivity of dust layers composed of individual 
dust aggregates 

Dust layers covering the surfaces of planets such as Mars or 
covering the surfaces of small bodies like comets or asteroids 
may not only be composed of single grains. It was recently dis- 
cussed by Skorov and Blum ( 201 1\ that the non volatile dust 
layer on the surface of a cometary nucleus could consist of in- 
dividual dust aggregates, which themselves are composed of 
micrometer-sized dust particles. These aggregates have formed 
before they were incorporated into the cometary nucleus. Fig. 
13 shows a sketch of a dust layer composed of dust aggregates. 
Here, the volume filling factor of the dust layer (piayer is a com- 
bination of the volume filling factor of the individual dust ag- 
gregates (pAgg and the volume filling factor of the packing struc- 
ture of the aggregates (l>s,ruc i- e. (piaye,- = (pAgg (psmic- 

In this Section, we show how to modify the model derived 
in Sect. 5.1 in order to calculate the heat conductivity of dust 
layers composed of individual dust aggregates. The total heat 




Contact Interface 
of the Monomere 



Contact Interface 
of the Aggregates 



R: Radius of the Aggregates 
f„: Radius of the Monomere 

Figure 13: Sketch of a dust layer composed of aggregates with radii R. The 
aggregates themselves consist of micrometer-sized dust particles with radii tq. 



conductivity of dust layers composed of monodisperse aggre- 



gates can be derived by modifying Eq. 12 and by taking the 



heat conductivity due to radiation into account (see Eq. |5]l, 

'^Layeriro, R, T, (pstruc, 4>Agg, MR)) - 



AAggifo, T, (pAgg) X 



9(1 -A-L) 



4£ 



^"'^ nyAgg(T)R^ 



1/3 



+ 'lRad(T,MR,(ps,ruc))- 



^(R, (pSiruc) 



(14) 



R, EAgg and jiAgg are the radius. Young's modulus and Poisson's 
ratio of the aggregates, respectively. The specific surface energy 
of the aggregates jAggiT) can be calculated, using the theory 
developed by |Weidhng et al.| ( |201 1] their Eqs. 5 - 7), 



jAggiT) = ^AggJ^Jio^iT) 



SiOi 



2/3 



(15) 



where ro is the radius of the monomeres. 



AAggiro, T, tpAgg) is the heat conductivity of the individual ag- 



gregates, which can be also calculated using Eq. 12 

9(1 



A-AggifO, T, (pAgg) = ^BulkiT) 



4E 



s/o. 



^ysioAnr^^ 



1/3 



X ^{rQ,4)Agg). 



(16) 



Note that the volume filling factor of the packing structure of 
the aggregates (psimc can be different to the volume filling factor 
of the individual aggregates (pAgg- 

For simplicity, we assume that the radiative heat transport 



(second term on the right hand side in Eq. 14 1 only takes place 
in the voids between the aggregates. Therefore, the heat trans- 
port inside the aggregates due to radiation is neglected. In this 
case, the mean free path of the photons can be calculated using 
Eq. 



13 and = ^st, 



Since dust layers and dust aggregates are not necessarily 
packed in one of the three close packing structures (see Tab. 
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Figure 14: Influence of the volume filling factor on the coefficient ^(r, (p) r (see 
Eq. [8]and|9}. The values for the simple cubic structure (sc; plus), the body- 
centered cubic structure (bcc; square) and the the face-centered cubic structure 
(fee; c ross ) are given by |Chan and Tien|jl973^ . A good fit to the data is given 
by Eq.fTslwith/i = (5.18 ± 3.45) X 10"^lnd/2 = 5.26 ± 0.94. 



|2]i, the coefficient ^(r, 0) (see Eq. |9]) lias to be modified in order 
to take otfier packing structures and, tfierefore, otlier volume 
filling factors of the material into account. The correlation be- 
tween the coefficient ^(r, 0) and the volume filling factor of the 
material was derived by fitting the dimensionless values (see 
Tab. [2]and symbols in Fig. 14 1 



^(r,0)r = 



1 



0.5315(0) Niir) 



r , 



with 



^(r,0)r = fi exp[/2<;*] 



(17) 



(18) 



14 1. The fit parameters are /i - (5.18 ± 



(dashed curve in Fig 

3.45)x 10^2 and /2 = 5.26±0.94. Extending the model with Eq. 
18 provides the possibility to calculate the heat conductivity of 
packed spheres for different volume filling factors (i. e. packing 
structures) of the material. 

The model for the heat conductivity of dust layers composed 
of individual dust aggregates was tested by a comparison with 
experimental results obtained by |Krause et al.| ( |2011) . In this 
work, the heat conductivity of dust layers consisting of individ- 
ual aggregates was investigated experimentally (their samples 
SI, S2 and S3, see Tab. [3]|. The dust layers were prepared by 
sieving the dust through metal filters with different mesh aper- 
tures. The mesh aperture determines the maximum size (diam- 
eter) of the aggregates. During these experiments, the surface 
temperature of the dust layer was between 300 K and 450 K. 



The experimental results are shown in Fig. 15 (symbols). Ad- 
ditionally, the most important properties of the samples and the 
used model parameters are summarized in Tab. |3] 

For the calculation of the heat conductivities of the sieved 
dust layers, we assume that the typical aggregate size is given 
by the mesh size. Furthermore, we used a Young's modulus of 
EAgg - 8 100 Pa, which was derived by Weidling et al. (2011 



Figure 15: Heat conductivity of dust layers composed of individual dust aggre- 
gates. The experimental results lI Krause et al.||20lT| symbols) are compared 
with the model for the heat conductivity of dust layers composed of individual 
aggregates (see Eq.Q4] curves). The calculations were conducted using the pa- 
rameters corresponding to the three different experiments (see Tab.[3j, 5 1 (blue 
dotted curves), S2 (red dashed curves) and S3 (green dashed-dotted curves). 
Two different temperatures were used for the calculations in order to take the 
lowest and the highest measured surface temperature of the dust layers into 
account, 300 K and 450 K. The straight lines are showing the calculated heat 
conductivity of the packed-aggregate structures without radiative heat transport. 
The model predictions, including radiative heat transport, ai'e visualized by the 
curves. Note that the experimental results can only be explained by the model, 
if radiative heat transport is taken into account. 



for aggregates with R - 0.5 mm. All aggregates, used in the 



work by Weidling et al. ( 201 l| l and in the work by Krause et al. 
( 201 l| l were composed of micrometer-sized Si02 particles with 
ro = 0.75 i^m. Thus, we also used this value for the calcula- 
tions. 

Fig. 15 shows the results of the experiments (symbols; 



[Krause et al.| [201 1| and the model (curves). The calculations 



were performed for two different temperatures, 300 K (low- 
est surface temperature measured during the experiments) and 
450 K (highest surface temperature measured during the experi- 
ments). For the calculation of the mean free path of the photons 
(see Eq. [TSj l, we used the volume filling factor of the packing 
structure of the aggregates. The variation of the aggregate sizes 
and the volume filling factors between the three samples leads 
to different coefficients ^(R, <p) R (see Eq. 18 and Tab. |3]l. 

Thus, we performed individual calculations for each exper- 
iment (see Fig. 15 1. For example, the blue dotted curves 



are showing the results of the model using the parameters of 
the Si sample. The red dashed and the green dashed-dotted 
curves correspond to the S2 and the S3 sample, respectively. 
The straight lines are visualizing the heat conductivity of the 
packed-aggregate structures without radiative transport. Ne- 
glecting the influence of the radiation, the heat conductivity de- 
creases with increasing radius of the aggregates. Obviously, the 
heat conductivity of the solid material (straight lines) can not 
explain the experimental results. 

However, taking the radiative heat transport into account 
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Table 3: Properties and model parameters for the dust layers used in the experiments conducted by 'Ki'ause et al.'l'2011V The maximum aggregate size is given by 
the mesh aperture (MA). Additionally, the volume filling factors of the aggregates <l>Agg, the volume filling factors packing structures (psimc, the corresponding mean 
free path of the photons inside the voids of the material A(i?) (see Eq. |13| and the coefficients f (ro, <f>Agg) ■ r and ^(R, >Psinic ) ■ R (see Eq. [Sjand Eq.|9j are shown. 



Experiment 


MA [mm] 


<l>Agg 


ipSmic 


HR,<l>Sln,c) [m] 


^(r(i,^Agg)- r 


i(R, <l>S,mc)-R 


Model 


Si (Plus) 


0.50 


0.45 


0.64 


0.8 i? 


0.55 


1.50 


Blue Dotted Curves 


St (Square) 


0.25 


0.45 


0.53 


1.2 i? 


0.55 


0.84 


Red Dashed Curves 


S3 (Cross) 


0.15 


0.45 


0.35 


2.5 i? 


0.55 


0.33 


Green Dashed-Dotted Curves 



(curves) leads to an increase of the heat conductivity for big- 
ger aggregates. If the radiative heat transport is dominant, the 
heat conductivity strongly increases with the radius of the ag- 
gregates. Hence, a comparison of our calculations with the ex- 
perimental data shows that the model can reproduce the exper- 
imental results within the uncertainty of the experiments (i. e. 
unknown mean temperature of the dust layers and experimental 
uncertainties, shown by the error bars.), if radiative heat trans- 
port is considered. However, the measured heat conductivities 
are showing a systematic trend relative to the model. The result 
of the S2 experiment fits quite well to the model (red dashed 
curves), while the result of the Si experiment can only be ex- 
plained with the calculation for 450 K (upper blue dotted curve) 
and the result of the S3 experiment is only in a good agreement 
with the model for 300 K (green dashed-dotted curves). This 
can qualitatively be explained by an increasing deviation from a 
linear temperature gradient from sample S3 through sample SI, 
due to an increased contribution of radiative heat transport (see 



also Sect. 6.2 and Fig. 18 1. The average temperature through- 



out the dust-aggregate layer increases from sample S3 through 
sample S 1 even if the bottom and surface temperatures are kept 



constant, which formally yields the tendency seen in Fig. 15 



It is important to note that the model has no free parameter 
to fit the theory to the experiments. The invented model only 
depends on the physical properties of the material, e. g. the 
grain size (pore size), the volume filling factor of the packing 
structure, the volume filling factor of the individual aggregates 
and the temperature of the dust layer. 

Furthermore, the experimental results can only be explained 
with the model, if radiative heat transport is considered, which 
consequently implies that radiation plays an important role for 
the heat transport in porous materials, even for temperatures 
T - 300 K. As a result, the heat conductivity of a porous ma- 
terial can not arbitrarily be reduced just by increasing the size of 
the aggregates or by just increasing the size of the monomeres, 
because of the influence of the radiative heat transport. 



0.010 



o 

CO 



tr 

CD 
X 



0.001 




1 



10 100 
Temperature [K] 



1000 



Figure 16: Temperature dependence of the Hertz factor (see Eq. |20( for parti- 
cles with a radius of r = 1 pm and a volume filling factor of the material of 
<4 = 0.5. 



6.1. Influence on the Hertz factor 

Here, we use the simple model (see Sect. |5.1| l to estimate the 
influence of the temperature of the material and of the radius of 
the particles on the Hertz factor H{r, T). The Hertz factor, 



H(r, T) 



(19) 



describes the reduction of the effective cross section of a porous 
material, due to the porosity of the material (see Sect. 4.3 1. 



Using Eq. 12 and 19 the Hertz factor is given by 



H(r, T) 



1 



4 E(T) 



^ ny{T)r^ 



1/3 



(20) 



6. Discussion 

In the first part of this section, we show how the simple model 
(see Sect. |5.1| l affects the Hertz factor. In the second part, we 
discuss the influence of the two models on the thermophysical 
modeling of cometary activity. Finally, in the last part of this 
section, we demonstrate the limitations of the models. 



Eq. 20 shows that the Hertz factor depends on the tempera- 
ture of the material and on the radius of the particles. 

Fig. [T6]shows the temperature dependence of the Hertz fac- 
tor for particles with a radius of r = 1 //m. Additionally, the 
influence of the grain size on the Hertz factor, for a tempera- 
ture of T - 300 K, is shown in Fig. 17 In both examples, the 
volume filling factor of the material is = 0.5. 
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Figure 17: Influence of the particle radius on the Hertz factor (see Eq. \20) 
for a temperature of 7" = 300 K and a volume filling factor of the material of 
<4 = 0.5. 





Figure 18: Sketch of a cometary surface layer. A solid H2O ice layer is covered 
by a porous dust layer. Additionally, the vertical temperature profiles inside 
the porous dust layer and inside the solid ice layer are shown, for dift'erent heat 
conductivities, A ^ MT) (porous dust layer; see Eq. |24| , A ~ (porous 
dust layer; see Eq. |26| and A ~ T"' (solid ice layer; see Eq. |28| . In this 
example, we assume that the heat conduction is dominated by radiation (A 4- 
A(J) <K /i ~ r^). Furthermore, the temperature of the cometary nucleus rA>„rf 
is fixed, for a better visualization of the influence of temperature dependence of 
the heat conductivity. Note that the surface temperature of the dust layer Tsmf 
and the temperature of the ice surface Tice are influenced by the dift'erent heat 
conductivities of the porous dust layer. 



6.2. Influence on the modeling of cometary activity 

The absorption of solar radiation by cometary surface layers 
leads to an increase of the surface temperature. Due to the heat 
conductivity of the material, heat can be transported into the 
interior of the cometary nucleus, which can cause sublimation 
of the icy constituents. While the plane albedo of the surface 
determines the amount of absorbed energy, the strength of the 
energy flux through cometary surface layers is given by the heat 
conductivity of the material and by the vertical temperature gra- 
dient inside the material. Modeling the energy distribution in 



cometary surface layers is a complicated task (see e. g. Davids- 
|son and Skorov| |2002| l, which is far beyond the scope of this 
work. However, we like to discuss the main implications of our 
model for the heat conductivity of dust layers on the thermo- 
physical modeling of cometary activity. 

In the simplest case, the structure of cometary near-surface 



regions can be divided into two layers. An ice-dominated layer 



and an ice-free dust layer, covering the ice (see Fig. 18 1. Then 



the surface temperatures of the dust layer and of the ice can be 
calculated using the energy balance equations at the surface of 
the dust layer and at the ice-dust interface. 



7(1 - A) = ecrr^V/ + ^D„.»(r, T, 0, A(r)) 



dT{K) 



dh 



(21) 



Surf 



and 



'^Dustir, T, 4>, A(r)) 



dT{h) 



dh 



Ice- Dust 



dT(h) 



dh 



+ Z(T„,)A(T„,). 



(22) 



ZiJice) is the sublimation rate of the ice and A(Tice) is the latent 
heat of the ice, which determines the cooUng of the ice due 
to the sublimation process. The vertical temperature gradients 
dT{h)/dh inside the dust layer and inside the ice can be obtained 
by solving the heat transfer equation. 

Computing the vertical temperature gradient for a constant, 
i. e. not temperature dependent heat conductivity, A + A(T), 
yields. 



dT(h) 



dh 



'Surf 



- T, 



Ice 



(23) 



^^-1(7") 



where, hjusi is the height of the dust layer. The corresponding 
vertical temperature profile inside the dust layer is then given 
by 



T(h)\ 



'Surf 



- T, 



Ice 



h +T,. 



(24) 



where h=0 is the position of the ice-dust interface. However, 
the assumption that the heat conductivity does not depend on 
temperature, is in contradiction with our model for the heat con- 
ductivity of dust layers (see Sect. [5] and Eq. 12 1. A tempera- 
ture dependent heat conductivity leads to a non-linear vertical 
temperature gradient inside the dust layer. Therewith, the total 
amount of transported heat through the dust layer is influenced. 
Thus, a constant heat conductivity should only be used, if the 
influence of the temperature on the heat conductivity is negli- 
gible, e. g. if radiative heat transport is not dominant or if the 
temperature differences across the dust layer are small (i. e. in 
< 
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our case for r < 15/im and T < 300 K; see Fig 
For the general calculation of the heat transport through porous 
dust layers, we propose to use the temperature dependent heat 
conductivity, introduced in Sect. [5] 

If the radiative heat transport is dominant, the temperature 
dependence of the heat conductivity is given by, A(T) ~ 
(see Eq. |5]l. In this case, the solution of the heat transfer equa- 
tion yields for the vertical temperature gradient inside the dust 
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layer, 
dT(h) 



dh 



1 



h 



X 



( '^Surf "^Ice ) 

-{TLf-TL) + Tl 



-3/4 



(25) 



Therewith, we can calculate the vertical temperature profile in- 
side the dust layer. 



T{h) 



A = ~ 



'Dust 



Surf 



- T 



Ice 



+ tL 



1/4 



(26) 



The heat conductivity of solid H2O ice also depend s on tern 
perature, AiceiT) - g\T'^, with = 567 Wm"' ( Klinger 



1980| l. Thus, the vertical temperature gradient and the vertical 
temperature profile inside the ice layer are given by 



dTih) 



dh 



7^1" — 

"Ice \ TnucI 



and 



X TnucI exp 



hice \ TnucI 



hice \ TnucI 



(27) 



(28) 



respectively. 

The influence of the different heat conductivities on the tem- 
perature profiles are shown in Fig. 
assumed A + A{T) «: A ~ 
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In this example, we 
(i. e. heat conduction is 
dominated by radiation). Thus, the amount of transported heat 
through the dust layer is increased from situation l(A AiT)) 
to situation 11 (A ~ T^). This leads to an increase of the ice 
temperature and to a decrease of surface temperature of the dust 
layer from situation I to situation 11. 

Tnuci denotes the temperature inside the nucleus, which is not 
known in detail. However, the internal temperature of cometary 
nuclei should be very low. One observational evidence is the 
presence of CO around the nucleus of Jupiter-Family Comets, 
which is interpreted as an indicator for very low temperatures, 
TnucI < 50 K ( |Thomas| |2009| l, because the subhmation tem- 
perature of CO ice is r w 26 K. Note that this interpretation 
can be misleading, because CO can also be trapped in the form 
clathrates inside other volatile components. 

It was shown that the usage of a temperature dependent heat 
conductivity influences the vertical temperature gradient inside 
the dust layer and, therewith, the amount of transported heat 
through the dust layer Thus, the usage of a temperature depen- 
dent heat conductivity for the calculation of the heat transport 
inside porous dust layers is crucial. 

Furthermore, the heat conductivity of a porous material can 
not be arbitrarily decreased by increasing the particle size, or 
by increasing the aggregate size (see Fig. [TTjand Fig. [TSj l. In- 
creasing the particle size, or aggregate size, automatically leads 
to an enhancement of the heat transport due to radiation. Con- 
sequently, the energy transfer on the ice surface can even be 
enhanced if the particle size is increased. 



6.3. Limitations of the model 

An underlying assumption of the model described in Sect. 
|5]is that the contacts between the particles or aggregates are 
formed by van der Waals forces. Thus, our model for the con- 
ductive heat transport in porous materials is limited to those re- 
gions in which the gravitational compression is negligible with 
respect to the adhesive force. The hydrostatic force per particle 
in the simple model is given by 



FHydr(r,h) = gpjrr h. 



(29) 



with g and p being the local gravitational acceleration and the 
mass density of the solid material, respectively. This compres- 
sive force must not exceed the van der Waals force (see Eq. [TO) l, 
if our model shall be applicable. 

Similarly, in the complex model, we get for the hydrostatic 
force per dust aggregate 

FHxdr,Agg(^^ <pAL, <Pdl) = g P h , (30) 

9 S true 

with 4'Surf - Psurf I P- Here, psurf is the mass density of the 



surface material. Following Wei dling et al.| ( |201 1 1, we get for 
the adhesive force per aggregate-aggregate contact 



Fv.dWAgg(R'T) = 3nyAgg(T)R, 



(31) 



where yAggiT) is given by Eq. [15] (see also Eq. 5 in Weidling 
et al. 201 1[ ). With suitable material parameters, we get y^gg = 
2 X lO^Jm-2. 

Using the above equations, we see that for the simple model 
and suitable g values for comets, F,.jii/(r, T) always exceeds 
FHydr{r,h) SO that our simple model (see Sect. 5.1 1 is widely 
useable. For example, Fydw{r,T) / FHydr{r,h) = 2 x 10^, for 
g - 3 X lO^'^ms"^, r = 1 mm, p = 1000kg m"-*, /i - Im 
andy = 0.02 J m^^. 



However, for the complex model (see Sect. 5.2 1, the hydro- 
static force exceeds the adhesive force for heights h > 0.5 m 
(for p = 1000kgm-^ ^ = 3 x 10""* ms^^, « = 1mm, 
ipsurf = 0.2 and (psirue - 0.5). Thus, our expressions for 
the conductive heat transfer in aggregates of dust aggregates is 
only valid for shallow depths. However, if the heat transport is 
dominated by radiative transfer, our model is always applicable. 

7. Conclusion 

In this paper, we investigated the heat transport in porous 
dust layers in vacuum. Therefore, the heat conductivity of Si02 
particles was measured experimentally (see Sect. |2]-|4|. Fur- 
thermore, we developed a model for the heat conductivity in 



loose granular media, based on the theory of Chan and Tien 
( [1973, see Sect.|5]l. 

This model was compared with the experimental results. We 
found that the model predictions are in a good agreement with 
the experimental results. Note that the model for the heat con- 
ductivity of porous dust layers has no free parameter to fit the 
obtained data. The model for the heat conductivity is capable to 
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distinguish between two different types of dust layers: dust lay- 
ers composed of single particles (simple model; see Sect. 5.1 



and dust layers consisting of individual aggregates (complex 
model; see Sect. |5.2| i. 

The heat conductivity of porous materials in vacuum is de- 
termined by two transport mechanisms, the heat conduction 
through the solid material and the heat transport due to radia- 
tion within the pores of the material. The heat transport through 
the solid material decreases with increasing particle size. The 
radiative heat transport strongly depends on the temperature 
of the dust layer and becomes important if the temperature of 
the dust layer is high (i. e. in our case for r > 15;um and 
T > 300 K; see Fig. [lO] and [TT]i. Additionally, both heat 



transport mechanisms depend on the packing structure of the 
material. 

In Sect. |6] we discussed the influence of the model on the 
Hertz factor. We found that the Hertz factor increases with tem- 
perature for fixed particle radius, due to an increase in the van 
der Waals strength. Furthermore, the Hertz factor decreases 
with increasing particle radius. Additionally, we discussed the 
implications of the model for the heat conductivity of porous 
dust layers on the thermophysical modeling of cometary activ- 
ity. We have shown that the temperature gradient inside the 
dust layer and, thus, the transported energy through the dust 
layer changes if the heat conductivity is temperature dependent 



6.3 



5.1 



and the com- 
if Fyi/wi'', T) exceeds 



(see Sect. 6.2 1. Thus, the exact determination of the properties 
of the dust layer, e. g. the particle size, the temperature and 
the volume filling factor, is important for the modeling of the 
energy fluxes inside cometary surface layers. 

Finally, we reviewed the limits of the model in Sect, 
Both models, the simple model (see Sect 
plex model (see Sect. |5.2| l, are useable 
FffrdrifJ'')- Using this condition, we found that the simple 
model is widely usable for suitable g values for comets. Fur- 
thermore, the complex model is also valid, if the height of the 
dust layer does not exceed h = 0.5 m. 
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